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We study the occurence of delay mechanisms other than 
periodic orbits in systems with time dependent potentials that 
exhibit chaotic scattering. By using as model system two har- 
monically oscillating disks on a plane, we have found the exis- 
tence of a mechanism not related to the periodic orbits of the 
system, that delays trajectories in the scattering region. This 
mechanism creates a fractal-like structure in the scattering 
functions and can possibly occur in several time-dependent 
scattering systems. 



Our study could serve as a first step towards the deeper 
understanding of certain scattering processes off oscillat- 
ing targets, such as molecules, atomic clusters, surfaces 
and nuclei in the cluster model. 

The outline of the paper is as follows: Tn Sec. we 
introduce our model. In Sec. [II we present the results of 
our numerical investigations. In Sec. IV we attempt to 



give an interpretation of the above results by identifying 
the mechanism that leads to the delay of the trajectories 
in the scattering region. In Sec. we summarize our 
main results and comment on the possible extensions of 
this work. 



I. INTRODUCTION 



II. DESCRIPTION OF THE MODEL 



Chaotic scattering has been an area of intensive re- 
search in recent years. The term chaotic scattering is 
used to indicate the very sensitive dependence of the vari- 
ables characterizing the outcoming trajectories on the 
variables characterizing the incoming trajectories [Q-Q. 
Chaotic scattering is usually associated with the appear- 
ance of fractal structure in the scattering functions. 

Although chaotic scattering has been extensively stud- 
ied in both classical and quantum mechanics, the vast 
majority of studies are concerned with systems where 
the scattering potential does not depend on time. A 
few studies of time-dependent chaotic scattering systems 
have recently appeared [|5| — |l~3|] but there is still a lack 
of understanding concerning the generic features of such 
systems. 

It is well known that in time-independent scattering 
systems the trajectories of the scattered particle are de- 
layed in the scattering region whenever the initial condi- 
tions of the scattered particle fall onto the stable man- 
ifold of one of the infinitely many unstable periodic or- 
bits of the system. This mechanism gives rise to sin- 
gularities in the delay time functions which are known 
to have a Cantor-set structure. The situation does not 
appear to be the same in time-dependent scattering sys- 
tems. Our aim is to investigate the occurence of delay 
mechanisms other than periodic orbits in time-dependent 
systems that exhibit chaotic scattering. We use a model 
system that is similar to that of Antillon et al. p2| . It 
consists of two circular disks on a plane that oscillate 
harmonically with time. The static counterpart of this 
system is integrable leading to smooth scattering func- 
tions 0. A harmonic oscillation law is used in order to 
avoid the superposition effects of the several harmonics. 



We consider the scattering of a free point particle off 
two circular harmonically oscillating disks on a plane, as 
shown in Fig. 1. The disks are considered to be much 
heavier than the scattered particle. The recoil of the 
disks at the collision is not taken into account and the 
energy of the system is not conserved. The centers of the 
disks are oscillating harmonically along the same axis. 
The position of the center of the i-th disk is given by the 
equation 



di(t) 



+ Aj sm(ujit + 4>i), 



(1) 



where denotes the equilibrium position of the center, 
A, is a vector directed along the axis of oscillation having 
magnitude equal to the amplitude of the oscillation, uji 
is the angular frequency of the oscillation and <f>i is the 
initial phase. The position r(i) and the velocity u(i) of 
the particle in the time interval between the (n — l)-th 
and the n-th bounce are given by 



r(t) = r„_i + (t - £„_i)u n _i 
u(t) = U„_i, 



(2) 



where r„_i denotes the position of the (n-l)-th bounce, 
u„_i denotes the velocity of the scattered particle after 
(n — l)-th bounce and denotes the time when 



the 



that bounce occured. In order to find the time t n when 
the next bounce will occur, we must solve the equation 



(3) 



where Ri denotes the radius of the i-th disk. This equa- 
tion must be solved twice (for i = 1, 2) and the smallest 
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non-negative solution has to be kept. The condition (|3|) 
leads to the equation 

ci(t - i n _i) 2 + c 2 sin 2 (cjji + <fo) 
+c 3 (t - t n -i) sin(u>it + (f>i) + Ci(t- tn-i) (4) 
+ c 5 sin(u!it + (pi) + c 6 = 

where the coefficients c± , . . . , cq are given by 

ci = u 2 _! 

c 3 = -2u„_i • Aj 
c 4 = 2(r„_ 1 -df ) )-u n _ 1 (5) 

C5 = 2(r n _ 1 -dj 0) ).A l 
c 6 = (r n _ 1 -df) 2 -^. 

Equation (4) is solved numerically to obtain t n . Using 
equation (|2|), the position of the next bounce is then given 
by 

r„ = r n _i + (t n - f„_i)u n _i (6) 

Since the mass of the disks is considered to be much 
larger than that of the scattered particle, the velocity of 
the particle after the n-th bounce is given by 

u„ = u„_x - 2[n ■ (u n _i - V»)] • n (7) 

where n is the normal to the disk at the point of the 
impact and Vj is the velocity of the involved disk at the 
instant of the collision, which is given by 

Vj = Aiuji cos(u!it n +<f>i). (8) 

In our study we have chosen the x-axis as the axis of 
oscillation (A, = (Ai,0)). We have also chosen the radii 
and the angular frequencies of the disks to be the same 
(i?i = R<2 = R, u)\ = uj2 = ui). By following the above 
procedure, we iterate numerically the map (r„, u n , t n ) — > 
(r„ + i, u„_|_i, t n+ i) and obtain the trajectories for a large 
number of initial conditions. 



III. NUMERICAL RESULTS 

A. Fractal-like scattering functions 

The scattering region is defined as a circular domain of 
radius Rq > Ri centered at the origin. As initial condi- 
tions we choose r(0) = and u(0) = (uq cos a, uq sin a). 
For the majority of our numerical calculations we have 
chosen the parameters of the system to be R = 1, 
d< 0) = -d< 0) = (-3/2,0), A 1 = A 2 = A = 0.1, 0! = 7T, 
4>2 = 0, lj = 10, a = I0~ 2 rad and Rq — 10. For these 
values of the parameters the scattered particle exhibits 
a small number of bounces (typically 4-7), in contrast 



to Ref. |12j where the number of bounces is much larger 
(typically more than 100) due to the different choise of 
the parameters. We have obtained the delay time T that 
the particle spends in the scattering region, the scattering 
angle (j> OM t and the outgoing velocity u out as a function 
of uq. The results are shown in Fig. 2. A rich fractal- 
like structure is observed in these scattering functions for 
rather small values of uq. As uq increases, the scattering 
functions become more regular. The dependence of T(uq) 
on the angular frequency to has also been studied. The 
results for u> = 100 and u> = 1000 are shown in Fig. 3. 
From this figure it is obvious that the change of u> simply 
introduces a scaling in the T(uq) function. In order to 
study the dependence of the structure of this function on 
the angle a, we have calculated T(uq) for a = 10~ 3 rad 
and a = lO^Vad. The results are shown in Fig. 4. We 
observe that as the angle increases the function becomes 
more regular in the sense that the range of uq in which 
T(uq) exhibits wild oscillations becomes smaller. In or- 
der to investigate the fractal structure of the system, we 
make successive magnifications of the T(uq) plot of Fig. 
2(a) in a region of uo around the value 0.8. The results 
are shown in Fig. 5. We observe that the fractal-like 
structure breaks at a very small scale. In this uo region, 
this breaking scale is found to be around 10~ 7 . In the 
following, a more in depth analysis of these findings will 
be given. 

B. Investigation of the fractal-like structure 

In order to quantify the observations concerning the 
fractality of the system made in the previous subsection, 
we have calculated the uncertainty dimension of T(uo) 
|l4],[l5]]. This quantity is given by d u = 1 — 0, where 
appears as /(e) « e 13 . /(e) is the fraction of uncertain 
points, for a given value of uncertainty, e, and for ran- 
domly chosen points uq. Each uq is considered to be un- 
certain if we find that the difference \T(uo) — T(uq + e)| 
is larger than a number of order 1. In our calculation 
we have included as many random points as necessary to 
obtain 150 uncertain points per run. We have calculated 
the uncertainty dimension for several values of a and u 
and in all cases log /(e) as a function of loge can be well 
fitted by a straight line. The results are shown in Table 
1. We can make the following observations: 

i) The dimension d u decreases with increasing angular 
frequency. This was expected since at the uj — > 
limit the system becomes regular. 

ii) The dimension d u decreases with increasing a. This 
was also expected since for smaller values of a the 
scattered particle senses more strongly the dynam- 
ics in the scattering region. 

For each angle a, there is a scale below which it is not 
possible to find any uncertain points. This supports our 
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observation that self-similarity does not persist at arbi- 
trarily small scales. 



C. Decay law 

We have also investigated the behaviour of the func- 
tion N(t)/No which gives the fraction of particles that 
remain in the scattering region after time t. It is known 
that for hyperbolic systems N(t)/No decays exponen- 
tially whereas for systems with marginal orbits and KAM 
tori it usually obeys a power law |16| . In order to stydy 
numerically the behaviour of N(t) /Nq we need to perform 
a Monte Carlo simulation by iterating a large number Nq 
of randomly selected initial conditions. In our system, in 
contrast to the static 2-disk system, an averaging over the 
initial velocity of the particle has also to be performed 
since the system is not conservative. In our calculation 
for the oscillating system, uq is uniformly distributed in 
(0, 2] since this is the range of uo where all the peaks of 
T(uo) occur. The angle a has been chosen in (0, 0.7] since 
for values of a in this range the particle exhibits at least 
one collision before exiting the scattering region. An av- 
eraging over the initial phases has also been performed. 
We have used 10 7 orbits. In order not to overemphasize a 
special class of peaks we choose Rq — 2.5 for our calcula- 
tions. The result is shown in Fig. 6. From this figure it is 
clear that for our system N(t) /N is very close to a power 
law with an exponent which for the chosen values of the 
parameters is found to be approximately equal to —2.38. 
The corresponding static system is expected to have an 
exponential decay law. In this case, no averaging in uq is 
performed since the system is conservative. If we average 
in mo for the static system we will also obtain a power 
law with an exponent equal to —2. If, on the other hand, 
we exclude from our calculations the "uncertain" trajec- 
tories (for which \T(uq,4>, a) — T(uq + e, (f>, a)\ is larger 
than a number of order 1) for the oscillating system we 
find a decay which for large times follows approximately 
a power law with an exponent equal to —5.09. The result 
is shown in Fig. 6. From the above it becomes clear that 
although the presence of the oscillation accelerates the 
escape of the particles, the presence of the high peaks in 
the T(uq) function (uncertain points) introduces a delay 
of the particles and slows down the escape. The origin of 
these peaks will be discussed in the following section. 



IV. INTERPRETATION OF THE RESULTS: A 
NEW DELAY MECHANISM 

By analyzing the trajectories that stay for long times in 
the scattering region, we found that these do not exhibit 
a large number of bounces. The observed delay comes 
from orbits along which the scattered particle loses much 
of its energy and therefore traverses segments of its orbit 
with a very low velocity. This sudden loss of energy can 



happen at any of the bounces provided that ||u„_i — 2Vi || 
is small but not as small as for the particle to bounce on 
the same disk again. We have classified the peaks of 
the T(uo) plot according to the bounce which leads to 
the major loss of energy. The result is shown in Fig. 
7. Prominent peaks that are due to the first bounce are 
observed. We also observe that around each of the peaks 
there is a rich fractal-like structure. In the following we 
are going to give a qualitative interpretation of the above 
observations by using a simple one-dimensional model 
which neglects the curvature of the disks. This model can 
yield quantitative results only for the first and the second 
bounce, since with the parameters chosen, the scattered 
particle senses strongly the curvature of the disks after 
the second bounce. 

Setting a = the dynamics is limited to one spatial 
dimension. In this case the system resembles to the Fermi 
acceleration model [[l7]jll|. The particle starts at the 
origin with initial velocity u(0) = (u ,0). In order to 
find the time t c when the first collision with the right 
disk occurs, we have to solve the equation 



Uot = — — R + Asm(u>t). 



(9) 



The time t c can be thought of as the abscissa of the 
first point of intersection of the straight line uot with the 
sinusoidal curve ^ — R + Asm(ujt) (see Fig. 8). From 
this figure it is obvious that t c is a discontinuous function 
of uq. Discontinuities occur for the values of uq for which 
the line is tangent to the sinusoidal curve. If we denote 
as t£ the value of t c after the discontinuity, the value of 
uq = Mq at which the discontinuity occurs is given by 



Uq = Auj cos(cjt+). 



(10) 



Therefore, the velocity of the particle after the first colli- 
sion u± is also discontinuous as a function of uq. At this 
point we should stress that the oscillation law need not 
to be only harmonic for these discontinuities to occur. 
In Fig. 9(a) a plot of u\ as a function of uq is shown. 
We observe sharp peaks that get denser as uq — * 0. For 
uq — > 0, the initial velocity that corresponds to a given 
t c is given by 



R — 



A 



Uq = 



tc 



(11) 



Therefore, the distance between two successive peaks on 
the uq axis is given by 



Au 



i Af 

tl 



(12) 



where At c = ^L. Combining the above two relations, we 
conclude that Auo ~ u\ as uq — > 0. The density of the 
peaks therefore increases as as uq — > 0. 

u a 

We observe that u\ gets close to near the discon- 
tinuities (uq = 2V, where V is the velocity of the disk 
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involved in the collision) . If the velocity after the colli- 
sion is small, but not small enough as for the particle to 
rebounce on the same disk, there is a delay of the particle 
between the first and the second bounce. We expect this 
delay to be present in the original system (a ^ 0) and 
to manifest itself as a peak in the T(uq) plot, since for 
small a the curvature of the disks can be neglected for 
the first collision. 

In the following we will give a qualitative interpretation 
of the fact that there is a rich fractal- like structure around 
the peaks of T(uq) for the original system . The quantity 
1 1 /«„ | is a measure of how much time is spent between 
the n-th and (n + l)-th collision for the one dimensional 
system. In Fig. 9(b) a plot of as a function of 

Mo for the one-dimensional system is shown. We observe 
that the peaks of this plot are very close to the peaks of 
the T(uq) plot that correspond to the delay of the particle 
between the first and the second bounce. The T(uq) plot 
can be thought of as several iterations of plots 
which are expected to have a structure similar to the 
\l/ui\ plot. For it that corresponds to a peak of the 
T(uq) plot, U\ falls in the low velocity region, where the 
structure of the |l/ui| plot is very dense. We therefore 
expect fingerprints of this dense structure to be apparent 
around any uo value which maps onto the low velocity 
region after some bounce. Since there is a lower bound 
on the velocity with which the particle can leave a disk, 
we do not expect the structure around the peaks to be 
infinitely dense. This is consistent with our observation 
that the fractal structure of T(uq) breaks at some small 
Mo scale. The higher a peak of the T(u ) plot is, the 
smaller the velocity of the particle along the orbit and 
the lower the breaking scale. The above mechanism also 
explains the presence of peaks in the scattering functions 
of Ref. |H| , however the structure there is more dense due 
to the different choise of the parameters. Furthermore, 
the mechanism is also present when the oscillating hard 
disks are replaced by two oscillating potential hills. A 
detailed analysis of the dynamics of such a system is left 
for a future study. 



V. CONCLUSIONS AND OUTLOOK 

In the present paper we have studied the effects of 
time-dependence on the scattering proccess of a point 
particle off two harmonically oscillating hard disks. A 
new mechanism leading to long-lived scattering trajecto- 
ries has been found. It is associated with the energy loss 
of the scattered particle at the collisions. This mech- 
anism is not related to the periodic orbits of the sys- 
tem and induces a fractal-like structure in the scattering 
functions. At the statistical level, the mechanism man- 
ifests itself as a change in the properties of the fraction 
N(t)/No of particles that remain in the scattering region 
after time t. Although this function still obeys a power 
law, the absolute value of the corresponding exponent is 



reduced by a significant amount. An interesting question 
that will be studied in the future is how time-dependence 
affects the dynamics of a system whose static counterpart 
is chaotic, such as the three disk system [^9|. As an ex- 
tension to this work, the transport properties of a lattice 
gas consisting of oscillating disks will be studied. An- 
other open question is the quantum manifestation of the 
new mechanism. 
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w(a.u.) 1 


10 


100 


a(rad) 








10" 3 


0.63 ±0.01 


0.62 ±0.01 


0.55 ±0.02 


10~ 2 


0.61 ± 0.02 


0.58 ±0.01 


0.54 ±0.01 


ht 1 


0.58 ±0.01 


0.53 ±0.02 


0.52 ±0.01 



TABLE I. The uncertainty dimension of the system for 
various values of the initial angle a and the angular frequency 
ui. 



FIG. 1. The two oscillating disks on the x — y plane. The 
positions of disks at time i n -i are drawn with the dotted line 
and at t„ with the solid line. The point A, which is defined 
with the vector r n _i, is the point of the (n — l)-th bounce. 
The point B, which is defined by r„, is the point of the n-th 
bounce. The dashed line with the arrows represents a segment 
of a trajectory that bounces between the two disks. 

FIG. 2. (a) The delay time T, (b) the scattering angle <j> ut 
and (c) the outgoing velocity a function of the initial 

velocity uo for angular frequency u) = 10 and initial angle 
a — W 2 rad. All the plots were created by iterating 10 4 
initial conditions. 



FIG. 3. The delay time T as a function of the initial veloc- 
ity u for q = 10" 2 rad and uj = 100 (a), u = 1000 (b). 



FIG. 4. The delay time T as a function of the initial ve- 
locity Uo for uj = 10 and a — 10~ 3 rad (a), a — 10~ 1 rad 
(b). 



FIG. 5. Magnifications of Fig. 2(a) around the value 0.8 of 
the initial velocity uo. Although in (a) a self similar structure 
appears to exist, it is found that this structure does not persist 
at arbitrarily small scales as shown in (b) and (c). 



FIG. 6. The fraction N(t)/No of particles remaining in the 
scattering region at time t. The thick solid curve includes 
all the trajectories and the thi n soli d curve includes only the 
"certain" trajectories (see Sec. Ill B). For large t both curves 
approximately obey power laws with different exponents. The 
dashed lines show the corresponding linear fits. The plots 
were created by iterating No = 10 7 initial conditions. 



FIG. 7. Some of the major peaks of Fig. 2(a) are classified 
according to the bounce that leads to the major energy loss 
of the scattered particle. The symbols A,D and o denote the 
first, the second and the third bounce respectively. 



FIG. 8. The abscissa of the point of intersection of the 
sinusoidal curve D/2—R+sin cut and the straight line uot is the 
time tc when the first bounce occurs for the one- dimensional 
system (a = 0). The slope of the straight line is Mo- For values 
of uo for which the straight line is tangent to the sinusoidal 
curve, t c is discontinuous as a function of uo- The value of t c 
after the discontinuity is denoted as tt ■ 
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FIG. 9. (a) The velocity ui after the first bounce as a func- 
tion of the initial velocity uo for the one- dimensional system 
(a = 0) . (b) 1 1 /tii | as a function of tto for the one- dimensional 
system. The locations of the singularities are very close to the 
locations of the peaks of Fig. 7 that correspond to a loss of 
energy of the scattered particle at the first bounce. 
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Abstract 

We study the occurence of delay mechanisms other than periodic orbits in 
systems with time dependent potentials that exhibit chaotic scattering. By 
using as model system two harmonically oscillating disks on a plane, we have 
found the existence of a mechanism not related to the periodic orbits of the 
system, that delays trajectories in the scattering region. This mechanism 
creates a fractal-like structure in the scattering functions and can possibly 
occur in several time-dependent scattering systems. 
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I. INTRODUCTION 



Chaotic scattering has been an area of intensive research in recent years. The term chaotic 
scattering is used to indicate the very sensitive dependence of the variables characterizing 
the outcoming trajectories on the variables characterizing the incoming trajectories |l]-f§|- 
Chaotic scattering is usually associated with the appearance of fractal structure in the 
scattering functions. 

Although chaotic scattering has been extensively studied in both classical and quantum 
mechanics, the vast majority of studies are concerned with systems where the scattering 
potential does not depend on time. A few studies of time-dependent chaotic scattering 
systems have recently appeared |5Hl3| but there is still a lack of understanding concerning 
the generic features of such systems. 

It is well known that in time-independent scattering systems the trajectories of the 
scattered particle are delayed in the scattering region whenever the initial conditions of 
the scattered particle fall onto the stable manifold of one of the infinitely many unstable 
periodic orbits of the system. This mechanism gives rise to singularities in the delay time 
functions which are known to have a Cantor-set structure. The situation does not appear to 
be the same in time-dependent scattering systems. Our aim is to investigate the occurence of 
delay mechanisms other than periodic orbits in time-dependent systems that exhibit chaotic 



scattering. We use a model system that is similar to that of Antillon et al. |12| . It consists of 
two circular disks on a plane that oscillate harmonically with time. The static counterpart of 
this system is integrable leading to smooth scattering functions [fj]]. A harmonic oscillation 
law is used in order to avoid the superposition effects of the several harmonics. Our study 
could serve as a first step towards the deeper understanding of certain scattering processes 
off oscillating targets, such as molecules, atomic clusters, surfaces and nuclei in the cluster 
model. 

The outline of the paper is as follows: In Sec. |I| we introduce our model. In Sec. 1 1 1 
we present the results of our numerical investigations. In Sec. [TV] we attempt to give an 



2 



interpretation of the above results by identifying the mechanism that leads to the delay 
of the trajectories in the scattering region. In Sec. [V] we summarize our main results and 
comment on the possible extensions of this work. 

II. DESCRIPTION OF THE MODEL 

We consider the scattering of a free point particle off two circular harmonically oscillating 
disks on a plane, as shown in Fig. 1. The disks are considered to be much heavier than the 
scattered particle. The recoil of the disks at the collision is not taken into account and the 
energy of the system is not conserved. The centers of the disks are oscillating harmonically 
along the same axis. The position of the center of the i-th disk is given by the equation 

di(t) = df ) + Ai sm{uit + (1) 

where df^ denotes the equilibrium position of the center, Aj is a vector directed along 
the axis of oscillation having magnitude equal to the amplitude of the oscillation, uji is the 
angular frequency of the oscillation and 0$ is the initial phase. The position r(t) and the 
velocity u(£) of the particle in the time interval between the (n — l)-th and the n-th bounce 
are given by 

r(t) = r n _i + (t- t n _i)u n _! 

u(t) = u n _ l} (2) 

where r n _i denotes the position of the (n-l)-th bounce, u n _i denotes the velocity of the 
scattered particle after the (n — l)-th bounce and t n -i denotes the time when that bounce 
occured. In order to find the time t n when the next bounce will occur, we must solve the 
equation 

- r(t)\\ 2 = Rl (3) 

where Ri denotes the radius of the i-th disk. This equation must be solved twice (for 
% = 1,2) and the smallest non-negative solution has to be kept. The condition ([5]) leads to 
the equation 

3 



d(t - t„_i) 2 + c 2 sin 2 ^ + <pi) 
+c 3 (t - £ n _i) sm(ujit + <pi) + c 4 (t - £ n _ x ) (4) 
+ c 5 sin^t + 0») + c 6 = 

where the coefficients ci, . . . ,cq are given by 

C l = U n-1 

c 2 = A 2 
c 3 = -2u n _i • Aj 
c 4 = 2(r n _i - dj 0) ) ■ u n _i (5) 

c 5 = 2(r n _ 1 -df ) )-A, 
c 6 = (r n _ 1 -df ) ) 2 - J R l 2 . 

Equation (4) is solved numerically to obtain t n . Using equation @, the position of the next 
bounce is then given by 

r n — r n~l + (t n — t n _l)u n _j (6) 

Since the mass of the disks is considered to be much larger than that of the scattered particle, 
the velocity of the particle after the n-th bounce is given by 

u n = u n _i - 2[n • (u n _i - Vj)] ■ n (7) 

where n is the normal to the disk at the point of the impact and V$ is the velocity of the 
involved disk at the instant of the collision, which is given by 

Vj = AiUi cos(uit n + (pi). (8) 

In our study we have chosen the x-axis as the axis of oscillation (Aj = (Ai,0)). We have 
also chosen the radii and the angular frequencies of the disks to be the same (Ri = R2 = 
R, u)\ =002 = 00). By following the above procedure, we iterate numerically the map 
(r n ,u n ,t n ) — > (r n+ i, u n+ %, t n+ i) and obtain the trajectories for a large number of initial 
conditions. 
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III. NUMERICAL RESULTS 



A. Fractal-like scattering functions 

The scattering region is denned as a circular domain of radius Ro > Ri centered at the 
origin. As initial conditions we choose r(0) = and u(0) = (wocosa, wosina:). F° r the 
majority of our numerical calculations we have chosen the parameters of the system to be 
R = 1, dS 0) = -d^ 0) = (-3/2, 0), A 1 = A 2 = A = 0.1, fa = vr, 2 = 0, u = 10, a = 10~W 
and Rq = 10. For these values of the parameters the scattered particle exhibits a small 
number of bounces (typically 4-7), in contrast to Ref. |12[] where the number of bounces 



is much larger (typically more than 100) due to the different choise of the parameters. 
We have obtained the delay time T that the particle spends in the scattering region, the 
scattering angle (p out and the outgoing velocity u out as a function of u . The results are 
shown in Fig. 2. A rich fractal-like structure is observed in these scattering functions for 
rather small values of uq. As uq increases, the scattering functions become more regular. 
The dependence of T(v,q) on the angular frequency uj has also been studied. The results for 
uj = 100 and uj = 1000 are shown in Fig. 3. From this figure it is obvious that the change 
of uj simply introduces a scaling in the T(uq) function. In order to study the dependence of 
the structure of this function on the angle a, we have calculated T(uo) for a = 10~ 3 rad and 
a = 10~ 1 rad. The results are shown in Fig. 4. We observe that as the angle increases the 
function becomes more regular in the sense that the range of u$ in which T{uq) exhibits wild 
oscillations becomes smaller. In order to investigate the fractal structure of the system, we 
make successive magnifications of the T{uq) plot of Fig. 2(a) in a region of uq around the 
value 0.8. The results are shown in Fig. 5. We observe that the fractal-like structure breaks 
at a very small scale. In this uq region, this breaking scale is found to be around 10~ 7 . In 
the following, a more in depth analysis of these findings will be given. 
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B. Investigation of the fractal-like structure 



In order to quantify the observations concerning the fractality of the system made in the 



previous subsection, we have calculated the uncertainty dimension of T(u ) fT^JT5|| . This 
quantity is given by d u = 1 — /3, where (3 appears as /(e) ~ e@. /(e) is the fraction of 
uncertain points, for a given value of uncertainty, e, and for randomly chosen points Uq. 
Each u is considered to be uncertain if we find that the difference |T(w ) — T{uq + e )| is 
larger than a number of order 1. In our calculation we have included as many random points 
as necessary to obtain 150 uncertain points per run. We have calculated the uncertainty 
dimension for several values of a and u and in all cases log /(e) as a function of log e can be 
well fitted by a straight line. The results are shown in Table 1. We can make the following 
observations: 

i) The dimension d u decreases with increasing angular frequency. This was expected 
since at the u — > limit the system becomes regular. 

ii) The dimension d u decreases with increasing a. This was also expected since for smaller 
values of a the scattered particle senses more strongly the dynamics in the scattering 
region. 

For each angle a, there is a scale below which it is not possible to find any uncertain points. 
This supports our observation that self-similarity does not persist at arbitrarily small scales. 



C. Decay law 

We have also investigated the behaviour of the function N(t)/No which gives the fraction 
of particles that remain in the scattering region after time t. It is known that for hyperbolic 
systems N(t)/No decays exponentially whereas for systems with marginal orbits and KAM 
tori it usually obeys a power law fl6| . In order to stydy numerically the behaviour of 
N(t)/No we need to perform a Monte Carlo simulation by iterating a large number No of 
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randomly selected initial conditions. In our system, in contrast to the static 2-disk system, 
an averaging over the initial velocity of the particle has also to be performed since the system 
is not conservative. In our calculation for the oscillating system, u is uniformly distributed 
in (0, 2] since this is the range of u where all the peaks of T(u ) occur. The angle a has 
been chosen in (0,0.7] since for values of a in this range the particle exhibits at least one 
collision before exiting the scattering region. An averaging over the initial phases has also 
been performed. We have used 10 7 orbits. In order not to overemphasize a special class 
of peaks we choose _R = 2.5 for our calculations. The result is shown in Fig. 6. From 
this figure it is clear that for our system N(t)/N Q is very close to a power law with an 
exponent which for the chosen values of the parameters is found to be approximately equal 
to —2.38. The corresponding static system is expected to have an exponential decay law. In 
this case, no averaging in u is performed since the system is conservative. If we average in 
uq for the static system we will also obtain a power law with an exponent equal to —2. If, 
on the other hand, we exclude from our calculations the "uncertain" trajectories (for which 
\T(u , 0, a) — T(w + e, 0, a) \ is larger than a number of order 1) for the oscillating system we 
find a decay which for large times follows approximately a power law with an exponent equal 
to —5.09. The result is shown in Fig. 6. From the above it becomes clear that although the 
presence of the oscillation accelerates the escape of the particles, the presence of the high 
peaks in the T(w ) function (uncertain points) introduces a delay of the particles and slows 
down the escape. The origin of these peaks will be discussed in the following section. 



IV. INTERPRETATION OF THE RESULTS: A NEW DELAY MECHANISM 

By analyzing the trajectories that stay for long times in the scattering region, we found 
that these do not exhibit a large number of bounces. The observed delay comes from orbits 
along which the scattered particle loses much of its energy and therefore traverses segments 
of its orbit with a very low velocity. This sudden loss of energy can happen at any of the 
bounces provided that ||u„_i — 2Vj|| is small but not as small as for the particle to bounce 
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on the same disk again. We have classified the peaks of the T(uq) plot according to the 
bounce which leads to the major loss of energy. The result is shown in Fig. 7. Prominent 
peaks that are due to the first bounce are observed. We also observe that around each of the 
peaks there is a rich fractal-like structure. In the following we are going to give a qualitative 
interpretation of the above observations by using a simple one- dimensional model which 
neglects the curvature of the disks. This model can yield quantitative results only for the 
first and the second bounce, since with the parameters chosen, the scattered particle senses 
strongly the curvature of the disks after the second bounce. 

Setting a = the dynamics is limited to one spatial dimension. In this case the system 
resembles to the Fermi acceleration model |17| , |18[ . The particle starts at the origin with 
initial velocity u(0) = (uq,0). In order to find the time t c when the first collision with the 
right disk occurs, we have to solve the equation 



The time t c can be thought of as the abscissa of the first point of intersection of the straight 
line uot with the sinusoidal curve ^ — R + Asm(ujt) (see Fig. 8). From this figure it is 
obvious that t c is a discontinuous function of u . Discontinuities occur for the values of u 
for which the line is tangent to the sinusoidal curve. If we denote as the value of t c after 
the discontinuity, the value of Uq = at which the discontinuity occurs is given by 



Therefore, the velocity of the particle after the first collision u\ is also discontinuous as a 
function of Uq. At this point we should stress that the oscillation law need not to be only 
harmonic for these discontinuities to occur. In Fig. 9(a) a plot of u\ as a function of uq is 
shown. We observe sharp peaks that get denser as uq — > 0. For uq — > 0, the initial velocity 
that corresponds to a given t c is given by 



u t = R + A sin(a>t). 



(9) 



Mq = Auj cos(u;£+). 



(10) 



u = 




(11) 



S 



Therefore, the distance between two successive peaks on the Uq axis is given by 

R TT 

Au = f 2 At c (12) 

where At c = ^j. Combining the above two relations, we conclude that Au ~ Mq as uo — ► 0. 
The density of the peaks therefore increases as as «o — > 0. 

We observe that u\ gets close to near the discontinuities (w = 2V, where V is the 
velocity of the disk involved in the collision). If the velocity after the collision is small, but 
not small enough as for the particle to rebounce on the same disk, there is a delay of the 
particle between the first and the second bounce. We expect this delay to be present in the 
original system (a ^ 0) and to manifest itself as a peak in the T(uq) plot, since for small a 
the curvature of the disks can be neglected for the first collision. 

In the following we will give a qualitative interpretation of the fact that there is a rich 
fractal-like structure around the peaks of T(uq) for the original system . The quantity 
is a measure of how much time is spent between the n-th and (n + l)-th collision for the one 
dimensional system. In Fig. 9(b) a plot of |l/wi| as a function of uq for the one-dimensional 
system is shown. We observe that the peaks of this plot are very close to the peaks of the 
T(uq) plot that correspond to the delay of the particle between the first and the second 
bounce. The T(uq) plot can be thought of as several iterations of |l/w n | plots which are 
expected to have a structure similar to the |l/tii| plot. For u that corresponds to a peak of 
the T(uq) plot, u\ falls in the low velocity region, where the structure of the |1/mi| plot is 
very dense. We therefore expect fingerprints of this dense structure to be apparent around 
any Uq value which maps onto the low velocity region after some bounce. Since there is a 
lower bound on the velocity with which the particle can leave a disk, we do not expect the 
structure around the peaks to be infinitely dense. This is consistent with our observation 
that the fractal structure of T{uq) breaks at some small Uq scale. The higher a peak of 
the T(uq) plot is, the smaller the velocity of the particle along the orbit and the lower the 
breaking scale. The above mechanism also explains the presence of peaks in the scattering 



functions of Ref. 12 , however the structure there is more dense due to the different choise 



of the parameters. Furthermore, the mechanism is also present when the oscillating hard 
disks are replaced by two oscillating potential hills. A detailed analysis of the dynamics of 
such a system is left for a future study. 



V. CONCLUSIONS AND OUTLOOK 

In the present paper we have studied the effects of time-dependence on the scattering 
proccess of a point particle off two harmonically oscillating hard disks. A new mechanism 
leading to long-lived scattering trajectories has been found. It is associated with the energy 
loss of the scattered particle at the collisions. This mechanism is not related to the periodic 
orbits of the system and induces a fractal-like structure in the scattering functions. At the 
statistical level, the mechanism manifests itself as a change in the properties of the fraction 
N(t) /No of particles that remain in the scattering region after time t. Although this function 
still obeys a power law, the absolute value of the corresponding exponent is reduced by a 
significant amount. An interesting question that will be studied in the future is how time- 
dependence affects the dynamics of a system whose static counterpart is chaotic, such as the 
three disk system 1 19fl . As an extension to this work, the transport properties of a lattice 



gas consisting of oscillating disks will be studied. Another open question is the quantum 
manifestation of the new mechanism. 
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TABLES 



ct(rad) 


cj(a.u.) 


1 


10 


100 


1(T 3 




0.63 ±0.01 


0.62 ±0.01 


0.55 ± 0.02 


io- 2 




0.61 ±0.02 


0.58 ±0.01 


0.54 ±0.01 


10- 1 




0.58 ±0.01 


0.53 ±0.02 


0.52 ±0.01 



TABLE I. The uncertainty dimension of the system for various values of the initial angle a and 
the angular frequency uj. 
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FIGURES 

FIG. 1. The two oscillating disks on the x — y plane. The positions of disks at time i n _i are 
drawn with the dotted line and at t n with the solid line. The point A, which is defined with the 
vector r n _i, is the point of the (n — l)-th bounce. The point B, which is defined by r n , is the 
point of the n-th bounce. The dashed line with the arrows represents a segment of a trajectory 
that bounces between the two disks. 

FIG. 2. (a) The delay time T, (b) the scattering angle (j) ou t and (c) the outgoing velocity u ou t 
as a function of the initial velocity uq for angular frequency uj = 10 and initial angle a = 10~ 2 rad. 
All the plots were created by iterating 10 4 initial conditions. 

FIG. 3. The delay time T as a function of the initial velocity uq for a = W~ 2 rad and u = 100 
(a), u = 1000 (b). 

FIG. 4. The delay time T as a function of the initial velocity uq for lu = 10 and a = W~ 3 rad 
(a), a = W~ 1 rad (b). 

FIG. 5. Magnifications of Fig. 2(a) around the value 0.8 of the initial velocity uo. Although 
in (a) a self similar structure appears to exist, it is found that this structure does not persist at 
arbitrarily small scales as shown in (b) and (c). 

FIG. 6. The fraction N(t)/No of particles remaining in the scattering region at time t. The 
thick solid curve includes all the trajectories and the thin solid curve includes only the "certain" 
trajectories (see Sec. Ill 1^ ). For large t both curves approximately obey power laws with different 



exponents. The dashed lines show the corresponding linear fits. The plots were created by iterating 
Nq = 10 initial conditions. 

FIG. 7. Some of the major peaks of Fig. 2(a) are classified according to the bounce that leads 
to the major energy loss of the scattered particle. The symbols A,D and o denote the first, the 
second and the third bounce respectively. 
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FIG. 8. The abscissa of the point of intersection of the sinusoidal curve D/2 — R+sincut and the 
straight line uot is the time t c when the first bounce occurs for the one-dimensional system (a = 0). 
The slope of the straight line is Uq. For values of uq for which the straight line is tangent to the 
sinusoidal curve, t c is discontinuous as a function of uq. The value of t c after the discontinuity is 
denoted as i+. 

FIG. 9. (a) The velocity u\ after the first bounce as a function of the initial velocity uq for the 
one-dimensional system (a = 0). (b) |l/iti| as a function of uq for the one-dimensional system. The 
locations of the singularities are very close to the locations of the peaks of Fig. 7 that correspond 
to a loss of energy of the scattered particle at the first bounce. 



15 



